Magnetic collective mode in underdoped cuprates: a phenomenological analysis 
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The dynamical spin susceptibility as relevant for underdoped cuprates is analysed within the 
memory-function (MeF) approach. A phenomenological damping function combined with a T- 
independent sum rule is used to describe the anomalous normal state and the resonant peak in the 
superconducting state, in particular its position and its relative intensity to the normal state. The 
relation with the random-phase approximation is discussed. The MeF method is generalized to the 
bilayer system in order to compare with inelastic neutron scattering experiments on YBa2Cu306+i 
which alows also for a quantitative comparison. In this context the problem of the missing integrated 
spectral intensity within the experimentally accessible energy window is also discussed. 

PACS numbers: 71.27.+a, 74.20.Mn, 74.25.Ha, 74.72.Bk 
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I. INTRODUCTION 

Since the first observation of the resonant mag- 
netic peak in superconducting (SC) optimally doped 
YBa2Cu306+a; (YBCO)fi the magnetic collective mode 
in cuprates and its role have been the subject of intensive 
experimental investigations, with the most direct infor- 
mation gained by the inelastic neutron scattering (INS). 
In YBCO it has been observed that with lower doping 
the resonant peak (RP) moves to lower energies, increas- 
ing at the same time in intensity^ In the last years, the 
intriguing feature of the RP hour-glass dispersion^ has 
been in the focus. Recently, also the existence of an opti- 
cal (even) resonant mode in the bilayer YBCO has been 
reconsidered pi It has been shown that the even mode can 
be well identified in underdoped (UD) regime, however 
with an essentially different intensity from the dominat- 
ing odd mode. On the other hand, it seems clear that the 
RP has to be closely related to the magnetic response in 
the normal (N) state, which has been well resolved by 
INS in UD cuprates. The response is typically that of 
an overdamped collective mode, however with an anoma- 
lous w/T scaling, being common both to UD YBCO^ and 
single-layer L^-zSr^CuO^i cuprates. 

Most theories, which address the RP and magnetic re- 
sponse in cuprates, are based on the treatment of the 
metallic system close to the antiferromagnetic (AFM) 
instability, describing the RP as a consequence of the 
AFM (over)damped soft mode in the N state and d x i_ y i 
SC gap in the electron-hole excitation spectrum, leading 
to a sharp RP below T c . Most frequently invoked is the 
RPA-like form for the dynamical susceptibility x q (<x>), de- 
rived in various ways or argued via the Hubbard, t-J or 
analogous models.- The latter ascribe the RP to a weak 
excitonic mode below the electron-hole continuum and 
seem to account qualitatively for optimum doping and for 
the overdoped regime, in particular for the position and 
the peculiar downward dispersion of the RP. 9 The alter- 
native approach, using the memory-function (MeF) de- 
scription, recently introduced by the present authorSfi2*ii 
gives analogous results in the latter cases, and in addition 



a more strongly pronounced upper dispersion 
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It is rather evident that the RPA-based theories for 
the magnetic response are less appropriate for the UD 
cuprates, even if taken in a broader sense as the phe- 
nomenological framework. In the first place, in the N 
state the usual RPA leads to T-independent and Fermi- 
liquid- like Xq(w), in contrast to the anomalous dynam- 
ics found by INS&i As discussed furtheron in more de- 
tail, the RPA-like description cannot account for a strong 
RP in UD cuprates and the spin-wave-like dispersion at 
higher energies, also observed in INS experimentspiiiii 
On the other hand, the MeF approach to X q (w) is able 
to give a unified description of the anomalous N-state 
as well as SC-state response, even more naturally in the 
UD regime. In this paper we will use the MeF formal- 
ism on the phenomenological level with some simplifying 
assumptions. It will be shown, that such an approach, 
generalized to a bilayer systems relevant for YBCO, can 
account for the T and doping evolution of Xq( w ) m the 
UD regime, in particular for the intensity of the coherent 
RP, and its jump at the onset of SC in relation with its 
position. The framework is also used to make a quanti- 
tative comparison with INS results for UD YBCO. 



The paper is organized as follows. In the following sec- 
tion we briefly sketch the MeF formalism (Sec. II. A) and 
then address the evolution of Xq( w ) with T, including the 
appearance of the RP and of its intensity relative to the 
sum rule C q (Sec. II. B). Next we present a critical com- 
parison of the MeF approach to the more commonly used 
RPA (Sec. II. C) and conclude the section with an exten- 
sion of the MeF approach to a bilayer system, as appro- 
priate to YBCO. In Section III we first use the available 
experimental data on YBCO6.5 and YBCOe.7 to extract 
the relevant parameters within MeF, and then present a 
quantitative comparison of the MeF approach to experi- 
ments. Conclusions are presented in Sec. IV. 
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II. MEMORY FUNCTION APPROACH 



A. Formalism 



generating always an undamped excitonic-like RP at 



LU r < UJ r . 



Within the memory function formalism the dynamical 
spin susceptibility can be generally written asiS 



ui 2 + wM q (w) — S q 



Ik 

Aq 



(1) 



where x q = x q (w = 0) and the 'spin stiffness' rjq — 

— i([-5f. q 'S'q)]) can be expressed with equal-time correla- 
tions for any microscopic model of interest and is in gen- 
eral a quantity only weakly dependent on q, T and even 
on hole doping c^>iS The latter is not the case for x q or 
for the 'mode' frequency w q . Instead of searching for an 
explicit approximation for either of these quantities, we 
fix them via the fluctuation-dissipation sum rule^S*ii 



^ Cth — X» = (S* q Sq> = C « 



(2) 



The underlying idea is that the equal-time spin corre- 
lations C q are much less sensitive on T. In particular, 
we conjecture that they do not change (significantly) at 
the N-SC transition, which we lateron verify for available 
YBCO data. Evidently, C q is q and doping dependent, 
but note that in strongly correlated systems C q is re- 
stricted by the total sum rule, (1/N) J^ q C q = (1 — c h )/A. 
In the following, we concentrate in the analysis to the 
commensurate Q = (tt, 7r), so the relevant quantity is 
C Q . 

Here, we do not intend to give a microscopic derivation 
for the damping function 7 q (w) = A/ q (cj), as performed 
within the t-J modeli 10 i 1:L i 12 Still, wc use the observation 
that the low-w damping within the N state of a doped 
AFM, being metallic and paramagnetic, is mainly due 
to the electron-hole excitations. If the Fermi surface of 
the doped system crosses the AFM zone boundary (as 
revealed by ARPES for most hole-doped cuprates) , in the 
N state one gets 7 q (w — > 0) > for q ~ Q. Hence, we use 
as the phenomenological input in the N state the simplest 
approximation 7 q (w) = 7, supported also by numerical 
calculations within the t-J models. The expression for 
X q (w), Eq.QJ, in the N state then reduces to a simple 
damped-oscillator form. 

In the SC state the d-wave gap introduces a gap also 
into 7 q (w). Here, we are not interested in the dispersion 
of the RPj2i±£ but rather in particular q = Q. Hence, we 
assume at T < T c 



7q(w<cj c ) = 0, 7 Q (cj>w c )=7, 



(3) 



where the effective gap uj c = 2A(q*) < A is given by 
the SC gap value at q* where the FS crosses the AFM 
zone boundary. For the discussion of the RP the corre- 
sponding real part is essential, 



M» = l]n|^| 

TT UJ C — UJ 



B. Qualitative analysis within the MeF 

The presented phenomenological theory has nontriv- 
ial consequences and predictions for the spin dynam- 
ics in cuprates, both for the N and thr SC state. It 
follows directly from Eq.Q) that the damped-oscillator 
form is well adapted to treat a collective mode in the 
UD regime close to the AFM instability. Here we focus 
on the q = Q mode and take into account also the fact, 
that in the N state the latter is generally overdamped, 
requiring 7 > ujq. Then, we get essentially two rather 
distinct regimes: a) if wq 3> ui c , consistent with modest 
Cq < 1, the N-state x'q^) IS ver Y broad leading to a 
weak excitonic-like RP at ui r < ui c . This is the situation 
corresponding to optimum-doped or overdoped cuprates. 

b) When luq is close or even below u) c , requiring Cq > 1, 
there is a pronounced low-w response already in the N 
state, transforming into a strong RP in the SC phase, 
exhausting a substantial part of the sum rule Cq . 

In Fig. 1 we present typical x'q( w )i corresponding 
to the UD regime. We fix the parameters from our 
knowledge of the results within the t-J model jiii as rel- 
evant for cuprates. For a weakly doped AFM we adopt 
77 ~ 2 J ~ Q.6t (note that t ~ 400 meV), the value for an 
undoped AFM. In the UD regime it is essential that the 
damping is small, 7 < t, both to get reasonable low-w 
response as well as underdamped spin waves at higher u>. 
Here, we assume 7 = 0.2i, while the SC parameters are 
chosen in accordance with data for UD cuprates, e.g., 
lu c = 0.15 t, T c = 0.05 t. Then C Q enters into x'qM 
as the only "free" parameter. In spite of simplifications, 
we observe in Fig. 1(a) several features consistent with 
experiments: a) the mode is overdamped in the N state 
with a Lorentzian form Xq( w ) lo/ (uj 2 + T 2 ) and a peak 
r oc T shifting with T, being the signature of the anoma- 
lous uj/T scaling, observed in UD cuprates&i b) Already 
at T > T c the peak T is below ui c , and in fact close to 
the position of the RP, which is not just a coincidence. 

c) At T < T c there appears a strong RP, exhausting a 
large part of the sum rule Cq. Also it is well shifted from 
the effective gap, i.e., Aw r = lo c — u) r is not small. The 
presented case is very close to the actual INS results for 
YBCO at ir = 0.5 - 0.7 j 2 i 5 i 13 as discussed lateron. 

We define the intensities as 



Jq = J g^x'qH- 



(5) 



(4) 



In Fig. 2 we follow the RP intensity I r , as a function 
of the relative RP position Au} r /uj c , where we vary Cq 
and fix the other parameters. In this way, we cover the 
span between the weak excitonic RP (small Cq) up to a 
pronounced AFM soft mode (large Cq). In INS experi- 
ments the strength of the RP is usually presented as the 
difference between the SC- and N-state response^ i.e., as 
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Figure 1: (Color online) Dynamical spin susceptibilities 
Xq(w) for different T above and below T c , with other pa- 
rameters fixed as given in the text, (a) Cq = 2, and (b) 
Cq = f . Spectra are additionally broadened with 5 = 0.01 t. 



Figure 2: (Color online) Resonant peak intensities I r (full 
lines) and AI r (dashed lines) , without and with N-state inten- 
sity subtracted, vs. relative position Au r /uc, for two values 
of 7 and varying Cq . 



AI r = I r (T ~ 0) - I r (T > T c 



(0) 



hence we plot in Fig. 2 also AI r , with subtracted N-state 
intensity at T > T c and Xq(^) integrated in the range 
< lo < u c . As discussed in Sec. IID, such a measure can 
be directly compared with the recent INS experiments: 4 
As noticed elsewhere^i& it follows from the singular be- 
havior of M'(Q, lj), Eq. (©, that I r oc Ac<j r , valid in the 
regime of weak RP, Aw r <C w c . From Fig. 2 one infers 
that a nearly linear dependence remains valid well be- 
yond this limit. Also the RP enhancement AI r shows 
similar dependence although overall reduced. Note that 
in the extreme case A7 r can become even slightly neg- 
ative. It is also worth observing that AI r vs. Auj r is 
only weakly dependent on 7, which might explain why 
YBCO with different doping x can be plotted roughly on 
a unique curved 



C. Comparison with RPA 

Let us also comment on the similarities and differences 
with the more usual RPA-like representation 



XqM = ftqM 1 - Jq] \ 

where in the SC state at T = 

1 \^f-t Ck+qEk + Ak+qAk 



(7) 



XqM 



^k+q^k 



Eh 



k+q 



^ - (£ k+q + £ k ) 2 ' 



(8) 



with Ek = v/^-F A^, whereas in the N state Xq cor- 
responds to the usual Lindhard function, and = 



—2t[cos(k x a) + cos(k y a)] — At' cos(k x a) cos(k y a) — fj, is 
the quasiparticle (QP) dispersion. The form, Eq.(J7Jl, has 
been derived in several ways or postulated based on mi- 
croscopic models such as the Hubbard model and the t-J 
modelfSsS, We treat it as a phenomenological expression 
where parameters jQ,t,t' are T independent, while still 
possibly doping dependent. 

In spite of the apparent different forms, RPA and MeF, 
Eq. Q, can be related at low ui. If the QP band crosses 
the AFM zone boundary, we get in the N state Xq(w — > 

0) = Tquj. Using the relation 

(xqM -1 - Jq)vq = S q -uj 2 - wMq(w), (9) 
it follows 



7 = Mq(^ 



0) = r/r Q /(x° Q ) 2 



(10) 



while r]Q = — [uj 2 xq]uj^oo oc t. Since in the SC state 
Xq(w < uj c ) = 0, the RPA result is also at T < T c 
formally close to our phenomenological MeF approach. 

Nevertheless, the qualitative and even more quanti- 
tative differences become very pronounced in the UD 
regime approaching the undoped AFM: a) in the N state 
Xq(w), Eq. |0, is essentially T independent and cannot 
account for the N-state anomalous scaling, b) For all 
reasonable QP bands, the effective 7 within the RPA ob- 
tained from Eq. (|10|) is very large, i.e., typically 7 > t. 
This prevents underdamped spin waves even at large w. 
c) The intensities I r of the RP are generally small, and 
in particular the N-SC difference AI r ~ 0. 

For illustration, we present in Fig. 3 characteristic re- 
sults within the RPA approach in the N and SC state 
T ~ 0, which correspond as close as possible to the 
regime of MeF results in Fig. 1. We choose the same 
u) c , while t = 3t' = 0.3i. Data in Fig. 3(a) correspond 
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Figure 3: (Color online)xq(<^) as evaluated within the RPA 
in the N state and SC state, respectively, (a) Jq = 0.95Jq, 
(b) Jq = 0.75J Q . 



to Jq = 0.95J q , i.e., close to the AFM instability. Al- 
though the RP is quite intense in this case, exhausting 
60% of the sum rule, the respective I r ~ 1.0 is con- 
siderably smaller in comparison to both the MeF case 
(Figs. 1(a), 2) and the INS data given later (see Table [IJ). 
On entering the SC state Cq exhibits a drop in mag- 
nitude, which may be quite substantial for Jq close to 
Jq, again in contrast with experiments. Similarly, for 
Jq = 0.75 J Q (Fig. 3(b)) where the RP intensity still ac- 
counts for ~ 20% of the sum rule, we get I r ~ 0.25, one 
order of magnitude smaller than the experimental data. 
Also note that AJ r remains negligeable (< 0.3) for all 

J Q < J Q- 



Table I: Values extracted from INS results for UD YBCO (odd 
mode). Cq and p xp (= 4//ig), are estimated from the data 
published in a) Ref.Q, b) Ref.El and c) Ref.U 





YBCO 6.5 


YBC0 6 .7 C 


T [K] 


5" 


85* 


100" 


12 


70 


200 


Cq[l/f.u.] 


2.1 


1.6 


1.9 


1.4 


1.4 


1.5 


/qV B /f.u.] 


26.4 


18.0 


11.3 


16.9 


14.9 


8.8 


/rV B /f-u.] 


19 b 






12.6 







where are related via the sum rules as before 



( Eq.© 

above interplane coupling JJl2, and we get 



Also, 77q'° can be explicitly expressed for the 



Ar; = r,° Q ~ Vq « Jx(8u ■ S 3i ) < [Jx/J)v- 



(13) 



The estimates in the literature are J±/J ~ 0.1 ji 5 - 
therefore to the lowest approximation we can neglect the 
difference, and t]q° ~ rj. The only appreciable source of 
distinction between the even and odd response is there- 
fore Cq > Cq. The latter difference can become sub- 
stantial in UD cuprates which are close to the AFM in- 
stability, where the soft mode is loq while Wq cx vJ±J 
remains gapped^ 

With the above simplifications we can apply previous 
single-plane results directly to bilayer YBCO. In partic- 
ular, in Fig. 2 results for I r ,AI r vs. Auj r for even and 
odd mode should fall on the same quasi-linear curve (for 
same 7, 77), corresponding to different Cq. E.g., we note 
that Aw r /w c = 0.05,0.3 correspond (for 7 = 0.2 t) to 
Cq = 1.2, 1.75, respectively, being qualitatively consis- 
tent with INS experiments. 5 



III. MeF APPROACH vs EXPERIMENTS 



D. Coupled- layers analysis 

In order to be able to describe more quantitatively the 
INS results for YBCO, which is a bilayer system, we gen- 
eralize the MeF approach to two coupled layers. We as- 
sume that the interaction between layers is only via an 
isotropic exchange^ 



H12 — J_L Sil • Si2- 



(11) 



Susceptibility can be defined as a matrix x l q as well 

as other quantities, 5^ ,r] 1 ^ etc. The plausible assump- 
tion we use here is that the dominant damping due to 
the electron-hole excitations can originate only from the 
intraplanar hopping, therefore M 12 = 0. Taking into ac- 



count the symmetry \ 



X etc. 



and odd functions, respectively, xtf 
can write in analogy with Eq. Q 



and defining even 
i\ 2 etc., we 



11 11 

Aq Ac 



lu 2 + ujMJw) 



re.o ) 



(12) 



Let us use the existing quantitative INS data for 
YBCO to check the relevance of the presented MeF 
analysis. First, the spectra for YBCO with x = 0.5i£ 
can be used to extract directly 7. Namely, at higher 
lo ~ 80 meV the spin-wave branches become under- 
damped and well resolved. From the width we can then 
estimate 7 ~ 70 meV ~ 0.2<, the value assumed in our 
presentation in Figs. 1(a), 2(a). We notice that such low 
7 (as well as spin waves) is essentially impossible within 
the RPA approach. 

Next, we extract data for the intensities Jq and the 
sum rule C Q , Eq. ©, for UD YBCO, both in the N 
and the SC state. To our knowledge such results are 
only available for x = 0.7— and x = 0.5. 2 Clearly, the 
spectra we can only integrate in the measured window 
lo < 70 meV, so results the represent a lower estimate. 

In Table U we present data for the odd mode, where 
the \aw-u) contribution is dominant. The important mes- 
sage is that Cq is indeed quite T-independent and in 
particular nearly conserved at the N-SC transition, be- 
ing the essential assumption within our MeF approach. 
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Figure 4: (Color online) MeF fit to INS data for XqM ( in 
absolute units) as given in Refs llsl Parameters are given in 
Table II. 



Table II: Values adopted in fitting INS data for UD YBCO 
(odd mode) as presented in Fig. 4. 





YBCO 6.5 


YBCO e.7 


T [K] 


5 


100 


12 


200 


rjQ [meV] 


130 


120 


130 


120 


7 [meV] 


35 


40 


60 


85 


cjq [meV] 


20 


26 


38 


40 


7 C [meV] 


30 




23 




C q [l/f.u.] 


2.1 


2.1 


1.2 


1.4 



Also, (theoretical) intensities Iq,I° are both large and 
comparable with MeF results in Fig. 2. 

Finally, we present an attempt of a quantitative fit 
of the same INS data2£ using our phenomenological de- 
scription as presented in Sec. IIA. We omit the YBCO6.7 
data at T = 85 K, which is just above T c as there seems to 
be already an indication for a reduced N-state damping, 
since the response appears already underdamped. In ad- 
justing the MeF form to the INS data we adopt the same 
form for Mq(w) as in Sec. IIA, but for T < T c we keep 
7q(o-> < uj c ) — 70 finite since experimentally the RP for 
T well below T c seems not to be resolution limited, i.e., 
the RP has a finite width. Thus: 



Essentially, r\ normalizes the intensities. While theo- 
retically (from model calculations) r\ t h ~ 2J is a very 
robust quantity, there are several reasons why rj t h be- 
comes renormalizcd. We are not dealing with the whole 
energy spectrum of spin fluctuations, but only with the 
low-frequency uj < J part. It is plausible, although the- 
oretically not well explored, that there is a substantial 
high-frequency uj > J dynamics, in particular increas- 
ing with doping. This is clearly inferred from spectra 
obtained in numerical simulations on finite clusters for 
the t-J model at finite temperature^ where at low T 
and low to moderate doping a rather distinctive sepa- 
ration of energy scales for uj < J and uj > It occurs. 
An enhancement of the sum rule Cq with respect to the 
Iow-oj region should then occur. But r\ should become 
enhanced even more since, being the first frequency mo- 
ment of Xq( w )j it is more susceptible to high-w tails in 
Xq( w )- Likewise, a partial reduction of the effective fluc- 
tuating spin would also result in a smaller r\. In any 
case, the same question already seems to be present when 
interpreting the INS results for undoped AFM where a 
reduced magnon intensity is observed through a renor- 
malization factor Z x ~ 0.5^2. 

Formally, a renormalization effect can be easily in- 
corporated into the MeF analysis by adding an effec- 
tive damping above some threshold uj* > J, assuming 
M£(w > uj*) ~ 7* > 7 and 



AMA(w) = — I In 

7T 



2A 7 
L 

TTUJ* 



(16) 



where A7 = 7* — 7. Such a term would renormalize 77 
for uj -c uj* 



77 — ► 77 = rj[l + 



2A 7l _i 



(17) 



which could explain the experimentally determined Z x . 
Note, however, that a separation of energy scales as men- 
tioned above would also lead to a similar result, since a 
suppression of Xq(w) in the energy interval uj* < uj < 2t is 
directly related to a local stepwise enhancement of damp- 



ing within the same energy interval, i.e., 
2t) - A 7 . 

IV. CONCLUSIONS 



M"[w* <uj < 



M'^uj) = 
M' q {uj) = 



1o0{uj c 



■hi 



uj) + j6(uj - uj c ) , (14) 

1=^1 • < 15 » 

UJr — UJ 



where 8(x) is the step function. In adjusting the MeF 
form we try to keep Cq,?7 = %j,7 as T-independent 
as possible. As before, we fix the effective SC gap at 
uj c ~ 60 meV, consistent with experimental analysis^. 
The resulting parameters are listed in Table II, while the 
fits are presented in Fig. 4. 

Overall, the fits appear quite satisfactory provided, 
however, that one chooses 77 ~ 120 — 130 meV ~ J. 



To summarize, we have presented a phenomenological 
analysis of the magnetic collective mode in UD cuprates 
within the MeF approach and have compared it with 
frequently used RPA in this context. Both approaches 
give some qualitatitvely similar results for the RP, in 
particular I r oc AclvA 1 ^ However, there are also clear 
differences. Due to large damping 7 ~ t, the RPA can- 
not capture the high-energy magnons^*i^ while these 
are easily reproduced within the MeF approach. Within 
RPA the RP intensities are too small in comparison to 
experiments, while Cq drops at the N-SC transition, in 
contrast with experiments as analysed in Table [I] 
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On the other hand, MeF shows in the N state an 
anomalous u/T scaling of Xq( w ); onr y being interrupted 
by the onset of SC, whereby the RP intensities I r close 
to experimental (see Table [IJ can be easily reproduced. 
However, there is significant discrepancy between AI r 
and AI° xp / A^i 2 B as the latter is quite small. In fact, in 
underdoped cuprates AI° xp might not be as relevant a 
quantity for it compares the RP (integrated) intensity be- 
low T c with the N-state intensity just above T c , where the 
appearance of a pseudogap can already reduce damping 
7 and induce a RP-likc responseAi^ii 

Generalizing the MeF approach to a bilayer system as 



appropriate for YBCO we have shown that the RP in- 
tensities in both channels, i.e., Ip° should fall on the 
same quasi-linear curve vs. Aw r , in agreement with 
experiments^ 

Finally, we have performed a quantitative fit within 
the MeF approach, which agrees well with the INS data, 
provided that r\ at low w is considerably renormalized. 

The authors acknowledge Y. Sidis and P. Bourges for 
the enlightening discussion. This work was supported by 
the Slovenian Research Agency under Grant P-106-017. 
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